Nuclear Spin Effect in Metallic Spin Valve 
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We study electronic transport through a ferromagnet normal-metal ferromagnet system and we 
investigate the effect of hyperfine interaction between electrons and nuclei in the normal-metal part. 
A switching of the magnetization directions of the ferromagnets causes nuclear spins to precess. 
We show that the effect of this precession on the current through the system is large enough to be 
observed in experiment. 



In recent years considerable theoretical and experimen- 
tal work is aimed at the controlled manipulation of elec- 
tron spin in nanoscale solid state systems, a field com- 
monly referred to as spintronics |]J. The main moti- 
vations for this research are applications in conventional 
computer hardware as well as the futuristic possibility 
of quantum computation 

3], 

using single electron spins as 
information carrying units (qubits). For both purposes, 
understanding the mechanisms of spin polarization, re- 
laxation, and dephasing in solid state systems is crucial. 

The branch of metallic spintronics has quickly evolved 
after the discovery of the giant magnetoresistance (GMR) 
in hybrid ferromagnetic normal metal structures |J, 
5]. Theoretical and experimental studies on magnetic 
mutlilayers have not only revealed interesting physics, 
but also already led to several applications in magneto- 
electronic devices. Magnetic recording read heads based 
on the GMR were first developed some ten years ago , 
but nowadays can be found in nearly all hard disk drives. 

In the context of quantum computation, semiconduc- 
tor quantum dots are regarded as promising candidates 
for storing electron spin based qubits Q ■ Recent progress 
in quantum manipulation of single spins |8( has over- 
come the effects of various spin relaxation processes in 
these devices. The unavoidable hyperfine interaction be- 
tween electron and nuclear spin presently attracts much 
attention. It has been identified as the main source of 
spin relaxation in high-purity samples at low tempera- 
tures 13, llOj and can even govern the electron transport 
in double dots At present, hyperfine interaction is 
seen as the main obstacle to demonstrate quantum com- 
putation with electron spins in solid state devices. 

In many other fields, for instance nuclear magnetic res- 
onance (NMR) experiments, hyperfine interactions play 
a central role already for decades. The Overhauser ef- 
fect is a common way to increase the degree of nu- 
clear polarization in metals enhancing NMR peaks. In 
semiconductors, optical orientation techniques |l.'j| are 
used to polarize the nuclear system. In the context of 
metallic devices, hyperfine interaction has been thought 
to be too weak to influence charge transport directly, and 
it has been regarded merely as an extra source of spin re- 
laxation 0. 

In this paper, we predict a clearly observable hyper- 



fine effect on electron transport in a metallic device. 
Thereby we demonstrate that hyperfine interactions may 
be important and possibly even dominant also for metal- 
lic spintronic devices. We consider electronic transport 
through a ferromagnet normal-metal ferromagnet multi- 
layer. This so-called spin valve is the basic magnetoelec- 
tronic device and the core component of all GMR based 
read heads. By changing the magnetization directions of 
the two ferromagnetic leads, one alters the total resis- 
tance of the device as well as the degree and direction 
of electronic polarization in the normal metal part in the 
presence of a current. 

Although the spin and particle transport properties 
of spin valves are well investigated and understood |l4| , 
effects of hyperfine interaction in magnetic multilayers 
have been hardly studied at all. One may think that 
these effects are negligible owing to the small value of 
the hyperfine interaction constant A ~ 10~ 6 eV in met- 
als. We show, however, that electron spins accumulating 
in the normal metal part can build up a significant po- 
larization of nuclear spins. The direction of this polar- 
ization is determined by the magnetizations of the leads. 
If the magnetizations are suddenly changed, this affects 
the electronic spin distribution in the normal metal part 
immediately (at a time scale r e ~ 1CP 11 s). The nuclear 
spin polarization reacts on a much longer time scale and 
will start to precess slowly around its new equilibrium 
direction. In this paper we are mainly interested in the 
feedback of the nuclear polarization on the electronic sys- 
tem. We show that due to such feedback the precession 
manifests itself as oscillations in the net current through 
the device. The amplitude of these oscillations is esti- 
mated as A/E t h- Here E t h is the Thouless energy char- 
acterizing the typical electron dwell time in the valve. 
The estimation is valid provided this time is shorter than 
the spin relaxation time r s f, which sets an upper bound 
for the effect, Ar s f /ft. For typical parameters, the rela- 
tive magnitude of the current oscillations can be of the 
order 10 -4 ~ 10~ 5 , which is clearly large enough to be 
measured in experiments. 

We model our system as a small metallic island con- 
nected to two ferromagnetic leads (Fig. We assume 
the island to be smaller than the spin diffusion length 
l s f and the time an electron spends in the island much 
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tors, defined as vE'(r) = Q~ l ' Y^k a ^ka e i where fl is 
the volume of the island. A n is the hyperfine coupling co- 
efficient between an electron and the nucleus at position 

f n , the vectors S n are the nuclear spin operators and a 
the Pauli spin matrices. 

We apply a second order perturbation expansion to 
find an expression for the time dependence of the elec- 
tronic distribution function f(k,t) = (fij. Q ,(£)afc j g(£)} 



N 



FIG. 1: A schematic picture of the system considered. A 
small metallic island is connected to two large ferromagnetic 
reservoirs with magnetizations m,L and rfiR. The contacts are 
charcterized by conductances Gl and Gr, which consist of 
spin-dependent parts G^ and G , and a mixing conductance 
G TJ \ The length of the normal metal part is significantly 
smaller than the spin diffusion length. 



smaller than t s f , which allows us to disregard spin-orbit 
relaxation mechanisms in the island. We also assume 
that the resistance of the junctions by far exceeds the 
resistance of the island itself. In this case, we can de- 
scribe the electronic states in the island with a single 
coordinate-independent distribution function f(E,t). 

The two ferromagnetic leads are modeled as large reser- 
voirs in local equilibrium, with magnetizations in arbi- 
trary directions and uir. Assuming for simplicity 
T = 0, we approximate the electronic distribution func- 
tion in the leads as fh(R){E) = #(a*l(,r) — E). The dif- 
ference in chemical potentials is due to the bias voltage 
applied eVj, = hl — We can disregard temperature 
provided eV& 3> fc^T. 

In our model, the electron spin polarization is mainly 
determined by the balance of spin-polarized currents 
flowing into and out of the ferromagnetic leads. However, 
a significant correction to this balance comes from hyper- 
fine coupling between the electron and nuclear spins. The 
resulting change of the polarization affects the net electric 
current in the device. So we will first derive an expres- 
sion for hyperfine induced polarization of electrons and 
nuclei, and then we combine the result with the known 
expressions for spin transport through spin valves. 

The Hamiltonian we use to describe the electronic and 
nuclear states in the island consists of an electronic part 
and a part describing the hyperfine interactions, 



H = H P 



hf 



H el = ^2 e k (4 T afcT + 4|«fcf) 

He 

H hf = 2 Aj n (t) ■ ¥(f n ) \<j *(f n ), 



(1) 



where are electron annihilation (creation) operators 
for spin up and down (a =?)!)■ We expressed the usual 
hyperfine contact Hamiltonian in electronic field opera- 



df(k,t) 
dt 



H(t),a{ a a kfj 
1 dt' 



H(t), H(t'),al a a k0 ), (2) 



where the indices a and j3 now span a 2 x 2 spin space. 
We see that the expansion can be completely expressed 

in the commuting operators S n and . 

Using Wick's theorem, we write the terms with four 
and six creation and annihilation operators as products 
of pairs, which then again can be interpreted as distribu- 
tion functions f(E,t). Further, we assume that the elec- 
trons are distributed homogeneously on the island and 
approximate A n — A/ no, A being the hyperfine coupling 
energy of the material and no the density of nuclei with 
non-zero spin . We find up to the second order 




A 



(S(t)) x f s (t) 



1 



- fs(t) ■ (S(t)) 



/.(*) 



(3) 



f P (t) [!-/*(*)] (S(t)), 



where we used the notation / = f, p t + f s ■ a, i.e. we 
split / in a particle and spin part. Nuclear spin enters 
here as the average polarization (S(t)). The first-order 
term describes the precession of electron spin in the field 
of the nuclei and disappears if electron and nuclear po- 
larizations are aligned. The second-order terms are all 
proportional to the hyperfine relaxation time defined as 
Thf — Tino/^A 2 v, v being the density of states at the 
Fermi energy. Since hyperfine interaction affects the elec- 
tron spin only, the contribution to the time derivative of 
f p is zero. 

This contribution is not the main one in the balance 
in the spin valve. Mainly, it is determined by electron 
transfers through the spin-active junctions. To describe 
this, we use the approach of 

Hi 

that is valid for non- 
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collinear magnetizations. This yields 



+ l{a L [1 - / p (t)] + /3 L m L • fat)} 
+ ^{ur \-f P {t)] + (3 R m R ■ fat)} 



rriL 
m R . 



(4) 



Following [14). we describe each spin-active junction with 
four conductances, G T , G x and G Ti = (G iT )*. (If the 
junction is not spin-active, G^ ~ G^ = G 1 ^ = G^ = G). 
The electron spin is subject to an effective field 



1 Im(G{ l )m L +Im(G R l )m R 



B„ 



Re(Gl i +G^) 



and we introduce dimensionless parameters characteriz- 
ing the spin activity of the junctions 



a L(R) 



L(R) 



G 



i 

L(R) 



2Re(G[ i + G R l 
2Re (GU„,) — G 



P, 



L(R) 



L(R) 



L(R) 



Lr L(fl) "T L(R) 



L{R) 



T _ r 4 

L( R ) ^L( R ) 



2Re(G[ x 



G R l ) 



The order of magnitude of (0J is estimated as l/r e , r e 
being a typical electron dwell time in the island, r e = 
e 2 M/2Re(G{ 1 +G#). 

Since by the essence of the spin valve, the spin balance 
affects the particle current through the device, the time 
derivative of f p is non-zero now, 



df p ^ 

dt ~T e { P L [ 



(a L m L + a R m R ) ■ f s (t)}. 



(5) 



We also need an equation for the time dependence of 
the nuclear spin polarization. One can obtain it from a 
perturbation expansion similar to J5J or directly inherit 
it from the fact that hyperfine interaction conserves the 
total spin of electrons and nuclei 



' d(S) \ 



dt 

/ hf 



-veV b 



dfs_ 
dt 



(G) 



hf 



We see from this that the relaxation time of nuclear spins 
Td ~ ThfUo/veVb- We assume that no other relaxation 
mechanism provides a shorther relaxation time. 

Let us now estimate and compare the time scales 
involved. For the nuclear system, the precession fre- 
quency uj = AveVb\fs\/finQ and the relaxation time 
Td ~ Tin^/irA 2 i> 2 eVb- Typical values for A range from 
10 _7 eV for weak coupling to 10 _4 eV (e.g. in GaAs (lif): 




we chose A = 5 x 10~ 6 eV. We take typical solid- 
state parameters to estimate no = 2.9 x 10 29 m -3 and 
v = 1.3 X 10 47 J _1 m -3 . For an applied bias voltage 
of Vb = 10 mV, this results in a precession frequency 
uj ~ 10 5 Hz and a nuclear relaxation time Td ~ 0.5 s. For 
the electronic system, the spin relaxation rate consists of 
two terms oc l/r e and oc I/t^/. We set the conductance 
of the F/N- interfaces to G ~ 3 SI -1 16] . If we choose 
dimensions of the metal island of 0.1 x 0.1 x 5 //m 3 , we 
find r e ~ 10~ n s. This corresponds to a Thouless energy 
Eth — fr/i~ e ~ 0.06 meV. The estimation for the hyperfine 
relaxation time reads Thf ~ 10~ 4 s. 

We conclude that on the time scale of all nuclear pro- 
cesses, f s and f p instantly adjust themselves to current 
values of voltage, magnetization, and importantly, nu- 
clear spin polarization. Their values are determined from 
the spin balance, 



= and ^=0. (7) 
dt 



As to nuclear polarization at constant voltage and mage- 
tization, it is of the order of 1 owing to a sort of Over- 
hauser effect produced by non-equilibrium electrons pass- 
ing the island. Indeed, it follows from Eq.^that the sta- 
tionary 2(5) = f s /[f p {l - f p ) + \f s \ 2 }. We see that (S) 
and f s are parallel under stationary conditions. This is 
disappointing since this will not result in any precession. 

The essential ingredient of our proposal is to change in 
time the magnetization(s) of the leads. Let us consider 
the effect of sudden change of the magnetization in one of 
the leads at t — 0. The electrons will find their new distri- 
bution, characterized by f new , on a timescale of r e . As we 
see from the nuclear spin system will start to precess 
around f new with the frequency estimated. The process- 
ing polarization will contribute to the effective field B ei 
B e —> B e + A(S(t))/h in Q. This will result in a small 
correction to f new , f^\ which is visible in the net cur- 
rent through the junction, due to its oscillating nature. 
A simple expression for this correction is obtained in the 
limit of weakly polarizing junctions («l,_r, Pl, R "C 1), 

f [1) {t) = ^{S{t))xf new oc A/E th , (8) 

while a more general expression is obtained by solving 
@ up to first order in A. The oscillatory part of the 
resulting current is given by 

I(t) = V b ^[P R m R - P L m L ] ■ /W(t). (9) 

The time dependence of nuclear polarization is still gov- 
erned by Eq. 

Combining Eqs |H1 and we find that the time- 
dependent current follows the behavior of (S(t)) and 
therefore exhibits oscillations with frequency u that are 
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FIG. 2: (a) Numerical calculation of the relative current fluc- 
tuations as a function of time. For t > 20 fis the i-axis is 
compressed. We used the parameters chosen in the text and 
further we took for both contacts a = 0.128, f3 — 0.115 and 
P = 0.333. The magnetizations switch at t — from m_t = 
(0, 0, 1), m R = (0, 0, -1) to m L = (0, 1, 0), m R = (0, 0,-1). 
(b) The dependence of the relaxation time and the frequency 
of the fluctuations on the asymmetry £ in the conductances, 
where f = hno / V A 2 'v 2 eVb and u> = AveVb/hriQ. 



damped at the long time scale t<j. The amplitude of these 
oscillations AI in the limit of small a and (3 reads 



\AI\_ 

(I) 



At 

-^■\(Soid) x /„ 



R m R 



- PL-mi 



(10) 



again proportional to A/Eth- In this equation (Sold) 
refers to the nuclear polarization before switching the 
magnetizations and fn^ R ^ is the part of mL(j?) perpen- 
dicular to fnew This relation makes it straightforward 
that one needs non-collinear magnetizations to observe 
any effect. 

In the same limit, the damping time and precession 
frequency are given by 



T d = T h f 



n 
veVb 



{e + r 1 + 2} 



(ii) 



ffTUfl of the amplitude is 4.4xlCT 5 . Eqs and U3) 
give Td = 0.63 s and lj = 1.0 X 10 5 Hz, in agreement 
with the plot. Typical currents through spin valves of 
these dimensions using a bias voltage of 10 mV range be- 
tween 10 and 100 mA. Oscillations of the order of 10 -5 
- 10 -4 should be clearly visible in experiment. The un- 
avoidable shot noise due to the discrete nature of the 
electrons crossing the junctions will not prevent even 
an accurate single-shot measurement, since the measure- 
ment time can be of the order of Td- An estimate using 
(5I) 2 ~ 2el/rd gives a relative error of 10~ 8 - 10~ 9 , at 
least three orders smaller than the oscillations. 

So far we have assumed precisely uniform electron dis- 
tributions. In a realistic situation however, the finite re- 
sistance of the island results in a voltage drop over the 
island, thus causing spatial variation of /o and f s . Impor- 
tantly, this gives variations in the precession frequency 
u> oc A\f s \/Ti. Such variation Aoj over the length of the 
island will contribute to an apparent relaxation of the 
spin polarization, since precession in different points of 
the island occurs with a slightly different frequency. This 
effect adds a term 1/r* = Aw to the damping rate l/r^. 
Assuming a simple linear voltage drop over the normal 
metal part, we find 1/r* — (Gj unc /Gj S i)ci;o, i-e. the ratio 
of the total conductance of the spin valve and the con- 
ductance of the metal island times the average oscillation 
frequency ujq. Although the effect can reduce the appar- 
ent relaxation time, provided (Acj)r^ 3> 1, it will not 
influence the time-dependent current just after t — 0. 

In conclusion, we have shown how hyperfine-induced 
nuclear precession in the normal metal part of a spin 
valve can be made experimentally visible. The precession 
should give a clear signature in the form of small oscil- 
lations in the net current through the valve after sudden 
change of the magnetizations of the leads. We found a 
coupled set of equations describing the nuclear and elec- 
tron spin dynamics resulting from a second order pertur- 
bation expansion in the hyperfine contact Hamiltonian. 
We presented a numerical solution for the net current and 
derived an estimate for the amplitude of the oscillations. 
We found that the relative amplitude of these oscillations 
is sufficiently big to be observable. 

The authors acknowledge financial support from FOM 
and useful discussions with G.E.W. Bauer. 



and 



AveV b \P R mji 



h n £ + 



(12) 



where £ = Gl/Gr characterizes the asymmetry of the 
conductance of the contacts. 

In Fig. [2k we plotted a numerical solution for the cur- 
rent I{t)/(I) and in[2) the dependence of id and Td on 
the asymmetry in conductance of the contacts. For[5R we 
made use of equations © and (7J , and inserted realistic 
&l.r and Pl.r- For the parameters used, the estimate 
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